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Abstract. We study the effect of dissipation on the infinite randomness fixed point 
and the Griffiths-McCoy singularities of random transverse Ising systems in chains, 
ladders and in two-dimensions. A strong disorder renormahzation group scheme 
is presented that allows the computation of the finite temperature behavior of the 
magnetic susceptibility and the spin specific heat. In the case of Ohmic dissipation the 
susceptibility displays a crossover from Griffiths-McCoy behavior (with a continuously 
varying dynamical exponent) to classical Curie behavior at some temperature T* . The 
specific heat displays GrifRths-McCoy singularities over the whole temperature range. 
For super-Ohmic dissipation we find an infinite randomness fixed point within the 
same universality class as the transverse Ising system without dissipation. In this case 
the phase diagram and the parameter dependence of the dynamical exponent in the 
GrifRths-McCoy phase can be determined analytically. 
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1. Introduction 

The interplay between quantum fluctuations and quenched disorder in the form of 
an extensive amount impurities or other random spatial inhomogeneities can lead to 
a new class of quantum phase transitions, governed by an infinite randomness fixed 
point (IRFP) as established for transverse Ising models [I] and many other disordered 
quantum systems (for an overview see [2]). Besides unusual scaling laws the transition 
the IRFP is characterized by a whole parameter range around the transition, in which 
physical observables display singular and even divergent behavior in spite of a finite 
spatial correlation length. This is the manifestation of Griffiths-McCoy singularities or 
quantum Griffiths behavior [HI [H [5l [6l [7]. They have their origin in rare regions of 
strongly coupled spins (or other quantum mechanical degrees of freedoms) that tend to 
order locally and thus produce a strong response to small external fields, long relaxation 
(or tunneling) times and small excitation energies. 

If the underlying quantum phase transition is governed by an infinite randomness 
fixed point the statistics of these rare events leads to a power law divergence of the 
susceptibility in a region around the quantum critical point with a continuously varying 
exponent. This dynamical exponent determines all singularities in the Griffiths-McCoy 
phase. Continuously varying exponents, interrelated in a specific way for different 
physical observables, were observed in many heavy-fermion materials, and it was argued 
that this is a manifestation of Griffiths-McCoy behavior due to an underlying IRFP 
[HIE]. In essence these systems form local moments that interact via long-range RKKY 
interaction and have a strong Ising anisotropy, such that an effective model describing 
these degrees of freedom and their interaction is a random transverse Ising system. 

Later it was argued that due to the interaction via band electrons the effective spin 
degrees of freedom are strongly coupled to a dissipative Ohmic bath [TOl[II]. From this 
point of view the rare regions should be described by spin-boson systems, which are 
known to behave classically for sufficiently strong coupling to the dissipative bath [T2] . 
which would destroy the expected Griffiths-McCoy singularities. 

Since in the presence of dissipation rare regions can undergo phase transitions and 
freeze independently from one another (like in the McCoy- Wu model in the mean-field 
approximation [I3]), the global phase transition of the system is destroyed by smearing 
because different spatial parts of the system order at different values of the control 
parameter [T4l [T5] . 

Recently we analyzed the random transverse Ising chain coupled to a Ohmic 
dissipative bath with a strong disorder renormalization group (SDRG) scheme and could 
demonstrate that the transition is indeed smeared, but argued that Griffiths-McCoy 
singularities are still observable, at least down to very low temperatures also in the 
presence of dissipation [H]. This was done by analyzing the gap and cluster distribution. 
In this paper we continue and extend this SDRG study by a) analyzing the low 
temperature behavior of the magnetic susceptibility and the spin specific heat in the case 
of Ohmic dissipation, where we will argue that Griffiths-McCoy singularities are visible 
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at all temperatures in the specific heat and above a (small) crossover temperature in 
the susceptibility; b) considering in addition to chains also ladders and two-dimensional 
systems, where we obtain similar results as for the chain; and c) applying the SDRG 
also to super-Ohmic dissipation, where we find a quantum phase transition belonging 
to the same IRFP universality class as the system without dissipation and compute 
analytically the phase diagram and dynamical exponent in the Griffiths-McCoy phase. 

The system that we study is the random transverse Ising model where each spin 
is coupled to a dissipative bath of harmonic oscillators, i.e. ferromagnetically coupled 
spin-boson systems [12]. It is defined on d-dimensional square lattice of linear size L 
with periodic boundary conditions (pbc) and described by the Hamiltonian 

(1) 



i ^ k 

where erf are Pauli matrices and the masses of the oscillators are set to one. The 
quenched random bonds Jj (respectively random transverse field hi) are uniformly 
distributed between and Jo (respectively between and ho). The properties of the 
bath are specified by its spectral function JTi{uj) 

J^i^) = ^ E —^(^ - ^k,) = la,n]-'u'9{Q, - u) , (2) 
^ k ^ 

where fii is a cutoff frequency and B{x) the Heaviside function such that 6'(x) = 1 if x > 

and 6'(a;) = if X < 0. The case s = 1 is known as Ohmic dissipation although s > 1 

(respectively s < 1) corresponds to a super-ohmic (respectively sub-ohmic) dissipation. 

Initially the spin-bath couplings and cut-off frequencies are site-independent, i.e. ai = a 

and = but both become site-dependent under renormalization. 

2. Real space renormalization. 

2.1. Decimation procedure. 

In this section, we derive in detail the real space renormalization scheme to study 
dissipative random transverse Ising model as in Eq. ([T]). For simplicity, we present the 
calculation in dimension d = 1 (extensions to higher dimensions are discussed below) 
and focus on the random transverse Ising chain (RTFIC): 



(3) 

k ^ ^ ^ 'i 

To characterize the ground state properties of this system (P), we follow the idea of a 
real space renormalization group (RG) procedure introduced in Ref. [17] and pushed 
further in the context of the RTFIC without dissipation in Ref. [3]. The strategy is to 
find the largest coupling in the chain, either a transverse field or a bond, compute the 
ground state of the associated part of the Hamiltonian and treat the remaining couplings 
in perturbation theory. The bath degrees of freedom are dealt with in the spirit of the 
"adiabatic renormalization" introduced in the context of the (single) spin-boson (SB) 
model [12], where it describes accurately its critical behavior [T8] . 



Finite temperature behavior of disordered quantum magnets with dissipation 



4 



2.1.1. When the largest coupling is a bond. Suppose that the largest coupling in the 
chain is a bond, say J2. The associated part H2 of the full Hamiltonian Hid in Eq. ([3]) is 

H2 



Hi'^ + V 



"2 "2 
T z z , f " z I Pk,i I 2 '^k,i 

-72^2^3 + 2^[Ck,^Xk,i(Ti + + i^k,i^ 



(4) 



V 



-h2al 



4=2,3 k 
x 

3 • 



Let us first focus on if2°'* in Eq. ([4]) and first introduce the notation for the spin part 
|S) = 15*2, 5*3), with Si = ±1 such that o^^\S) = Si\S) for i = 2,3. Considering now the 
two baths on site 2 = 2,3 respectively, they are composed of a set of harmonic oscillators 
which are labeled by an integer k (which formally runs from to 00) and by i = 2, 3. 
We denote by \nk,i) the eigenvalues of these harmonic oscillators such that 



Pk,i . 2 "^k,' 



\nk,i) = {nk,i + -)uJk,i\nk,i) 



(5) 



H2^^ can be 



In the absence of the coupling between the spins and the bath 
straightforwardly diagonalized by tensorial products of |S) and \nk,i). The corresponding 
eigenvalues are simply the sum of the eigenvalues of the individual Hamiltonians in Eq. 
dH) without the last term of interaction. The coupling between the spins and the baths 
does not change these eigenvalues (up to a global shift) and only affects the eigenstates. 
We introduce the shifted "eigenvectors" 



\n 



) = exp 




J 



2,3 



(6) 



from which we can construct the eigenvectors of H2 
the baths and the spins as 



(0) 



S,n) 



n. 



including the interaction between 



(7) 



The eigenvalues of H2^^ are given by 



H§'^\S,n)=EZ\S,n) 



(0) 



E, 



(0) 
S,n 



J2^2^3 + E E 

i=2,3 k 



{nk,i + 



1 ^kA 



2 



(8) 



Each level is thus a priori degenerated twice (except accidental degeneracy) and in the 
limit of large coupling J2 we first restrict ourselves to the lowest energy levels, such 
that S'2S'3 = +1. Performing perturbation theory in V , one obtains that the first order 
corrections vanish. To second order in V , one has to diagonalize the 2x2 matrix V^*^^) in 
the eigensubspace associated to the zeroth order eigenvalue with 5*2 5*3 = +1 which 
is formally given by 

r|S',n')(S',n'|V^ 



y{2) 



E 



E. 



(0) 
S,n 



E. 



(0) 
S',n' 



(9) 
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One obtains from ([9]) the diagonal elements 



In'- 



2J2 + Efc(^fc,2 - nU)c^fc,2 

"2 ' 



and the off-diagonal elements 

VlV = V21 = ^('^2 |t^2 )('^3 K) ■ (11) 

In the absence of a coupling to the dissipative bath {i.e. Ck,i = for all k and i) 
the shifted eigenstates ((61) are identical with the non-shifted eigenstates and therefore 
(nr|n+) = and thus V^^^ = v!^ = and vl^^ = l^?^ = -h2h^/ J2. This matrix 

has two eigenvalues whose difference, the gap, is 2/2,2^3/ J2- Thus for each oscillator state 
the low lying excitations of if2 in ((11), with 82 = S^, can again be described by an effective 
two-state system, i.e. a spin in a transverse field of strength h' = /i2^3/<^2- The spirit of 
the strong disorder renormalization group is to keep this effective two-level system (for 
each oscillator state) and to neglect the large energy doublet with S2 ^ S3. In this way 
one has replaced two spins (with moments fi2 and /is) and a large coupling J2 between 
them by a single effective spin with moment /i2 + /is in a small transverse field /i2^3/</2, 
thus one degree of freedom with a large energy has been decimated. 

In the presence of non-vanishing couplings Ck,i to the oscillators one needs to 
decimate also the high energy modes of the bath such that Uk^i > pJ2, where p is 
some (large) number. Given that J2 is a large energy scale, the low lying energy levels 
will be those with n^ j = for u;^ j > pJ2. Therefore we decompose the oscillator states 
according to 

|nf) = |n<^) ® |n>^) , (12) 
with S* = ±1 and 

\<') = <S>K) and |n>^) = (g)|nfj, (13) 

k< A;> 

where kf = {k\ujk^i < PJ2} and kf = {k\ujk^i > pJ2}- Additionally we introduce the 
product state of oscillators which are in the ground state by |of ') = 0;j> |0f^). At 
energy scales smaller than pJ2 all oscillators with frequencies larger than pJ2 will be in 
their ground states, and therefore we will consider the matrix elements in (fTOl) and (fTB 
only for oscillator states In"*") = |n'^"^) ® \o^^)- For these states the two sums on the 
r.h.s. of (fTOl) read 



To leading order in J2 one can neglect the term J2k<{'^k,2 — ^k,2)^k,2, since it involves 
only frequencies smaller than pJ2. Then the sum over the low frequency oscillator states 
|mf ) yields one since they form a complete basis for the low frequency oscillator Hilbert 
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space and the individual terms in the sum do not depend on the quantum numbers nf 
any more. Thus the diagonal matrix elements in (fTOl) read to leading order in J2 



v!P-r^ l(o?>r-)P l(o?>.?')|- , ,,,, 

Note that this expression does not depend on the quantum numbers n< for the low 
frequency oscillators. For the non-diagonal matrix elements in ( fTTll one gets 

V{? = Vi? = -^(n<-|n,-)(n3<-|n3<+) (16) 
■-'2 

with 

^= (o>+ior)(orio3>-) = n(o^jo,-)n(o+ io,-) . (i?) 

The amplitude A can be then expressed in terms of the spectral density, using that 
(Ofc,i|Ofc,i) = exp-iClJui;). This yields 

^ = exp / -Ar^du - - / -^du , 18 

\ 71 JpJ2 TT JpJa j 

where we have used the definition of the spectral density in Eq. ([2]). Since the diagonal 
term does not depend on n*- the diagonalization of V"^^^ yields (up to second order) the 
following correction to the lowest eigenvalues 

1=2,3 k ^ -^2 

(19) 

where = —J2 + Vn — | Z]i=2,3 Sfe Cli/i^fc « is a constant, independent of n<. We now 
consider an effective spin-boson Hamiltonian coupled to both baths 2 and 3: 

H2 = -hal + E E (4^ + + CkAk^^ri) (20) 

where the frequencies are such that ujk,i < pJ2- The effective spins being coupled to 
both baths, one has 

MUJ) = e{pJ2 - U){J2{UJ) + MUJ)) . (21) 

Treating the small parameter h2 in (degenerate) perturbation theory, one obtains the 
low lying eigenvalues of H2 to first order in h2: 

ESn< = E E i^k, + I) ujk, ± /^2(nr|n2<+)(n3<-|n3<+) . (22) 

1=2,3 k ^ 

The comparison between Eq. ffTOl) and Eq. ([22l) shows that the low energy spectrum of 
the two interacting spin-bosons in H2 can be described by a single spin-boson system 
with renormalized parameters given by 
r Ah2h^ 

h2 = — z — (23) 
'J2 

(52 = ^2 + "3 , ^2 = PJ2 , (24) 
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where A, which depends on the parameters of the Hamiltonian H2, is given by Eq. (fTSll 
and the equalities in (l24l) are a direct consequence of Eq. (l2Tl) . This effective spin-boson 
will interact ferromagnetically with the spin-boson on site 1 and site 4 with couplings 

Ji = Ji , J2 = J3. (25) 
These relations ( l23l ). (l24l) and (l25l l constitute the first set of decimation rules. 

2.1.2. When the largest coupling is a field. Suppose now that the largest coupling in 
the chain is a transverse field, say h2. Before we treat the coupling of site 2 to the rest 
of the system — Jicrfcxl — J2(Tlal perturbatively as in [3] we consider the part of the 
Hamiltonian that represents a single spin-boson system 

H'^ = -h2al + ^(^ + a;,,2^ + Ck,2ik,2(Tl) • (26) 

k 

In this case, one would like to have a way to decimate the high energy modes of the 
bath, here the harmonic oscillators such that 0;^ 2 > P^2, where p is some (large) number. 
Since for those oscillators cj^ 2 ^ ^2 one can assume that they adjust instantaneously to 
the current value of erf the renormalized energy splitting is easily calculated using Eq. 
(I22II - the so called adiabatic renormalization [12]- and one gets an effective transverse 
field h2 < h2: 

h2 = A'h2 , ^2 = ph2 (27) 
A' = e.J-lf"'^,.). (28) 

y TT Jph2 ) 

If h2 is still the largest coupling in the chain the iteration f l27ll is repeated. Two situations 
may occur depending on the parameters s and a^. 

• Ifs>lors = l and 0:2 < 1 this procedure (l27ll will converge to a finite value h\ 
given by 



/in = /i2 exp / — T^dv ~ { 

\ TT J phi j 



/i2exp ( -j-^ ) ) s>\ 



' ph 



i2_ 



^2 \ ^1 , S = 1,02 < 1 



(29) 



where the expression of h*2 for s > 1 is valid only in the limit p/12 <^ ^^2- In this 
case, the spin-boson system at site 2 is in a delocalized phase in which the spin and 
the bath can be considered as being decoupled (formally 0:2 = 0), as demonstrated 
by an RG treatment in [H]. If this value h\ (1291) is still the largest coupling in 
the chain the spin on site 2 will be aligned with the transverse field. As in the 
RTFIC without dissipation, this spin is then decimated (as it will not contribute to 
the magnetic susceptibility) and gives rise, in second order degenerate perturbation 
theory, to an effective coupling J\ between the neighboring moments at site 1 and 
312] 
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• If s = 1 and 0:2 > 1, /i2 can be made arbitrarily small by repeating the procedure 
([27ll implying that the SB system on site 2 is in its localized phase [18] and 
essentially behaves classically: the decimation rule (l27ll indeed amounts to set 
/12 = 0. Such a moment, or cluster of spins, will be aligned with an infinitesimal 
external longitudinal field and is denoted as "frozen". 

These relations f l27ll30ll constitute the second set of decimation rules. The complete 
decimation procedure is sketched for the Ohmic case in Fig. 1. 




Figure 1. Left: Sketch of the adiabatic renormahzation of the bosonic bath for 
a single spin-boson system in the Ohmic case: The fast oscillators with frequencies 
uj > ph, where h is the transverse field acting on the spin and p an arbitrary large 
parameter, are treated in the adiabatic approximation resulting in a renormalized 
transverse field h for the spin and a decreased upper cut-off frequency f2max for the 
bosonic bath. This procedure is continued until it runs into a fixed point, where either 
the renormalized field h* vanishes (for a > 1) and the spin is frozen, or the renormalized 
field h* has a positive value, in which case it is a nonfrozen spin in a transverse field h*, 
decoupled from the bath. Right: Sketch of the SDRG in the presence of a dissipative 
bath for the Ohmic case. When a field h is the largest coupHng (top), first the oscillator 
bath is renormalized along the lines sketch on the left site. Once this is done and the 
field is still the largest coupHng the spin gets decimated as described in the text. When 
a bond J between two spins at site 2 and 3 is the largest coupling the two spins get 
decimated to form a cluster with moment fi2 + /^3, coupled to a combination of bath 
2 and 3, i.e. a new dissipative bath with coupling strength a2 + a^. 



2.2. Numerical implementation. 

In the following we analyze this RG procedure defined by the decimation rules (I23II251 \27\ - 
[30I I numerically. This is done by considering a finite system of linear size L with pbc and 
iterating the decimation rules until only one site is left. This numerical implementation 
has been widely used in previous works [2], [19], and it has been shown to reproduce 
with good accuracy the exact results of Ref. [3] for the RTFIC [Uj. In particular, the 
transverse field h acting on the last remaining spin is, at low ferromagnetic coupling Jo, 
an estimate for the smallest excitation energy. Its distribution, Pl^H/Tq), where Fq is 
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the largest coupling in the initial system of linear size L, reflects the characteristics of 
the gap distribution [5]. This quantity, and specifically its dependence on the system 
size L can be efficiently used to characterize Griffith-Mc Coy singularities and critical 
behavior characterized by an infinite randomness fixed point. 

3. Ohmic dissipation. 

Ohmic dissipation means s = 1 in ([2]), «. e. a spectral function for the oscillators that is 
linear in the frequency (up to the upper cut-off For a single spin in a transverse field 
h and coupled to such an Ohmic bath, a lot of results are available [I2]. Here we mention 
only that this system has a phase transition at zero temperature driven by the coupling 
strength a. For small a the spin can still tunnel quantum mechanically, whereas for 
large a the spin is frozen and behaves classically, the critical coupling strength adh) 
is equal to 1 in the limit where h/fl 1 where fl is the cut-off frequency of the bath, 
an exact result predicted correctly by the adiabatic approximation mentioned above. 
Such a transition is also present in an infinite ferromagnetic spin chain coupled to a 
dissipative bath, as it was shown recently numerically [201. Here we want to focus 
on the interplay of disorder, quantum fiuctuations and dissipation and study random 
transverse field Ising systems coupled to a dissipative environment by implementing the 
decimation rules (12311251 [27ll30ll for ohmic dissipation. For s = 1, the amplitudes A in 
Eq. (fTSl) and A' in Eq. (l28ll which enter these decimation rules are given by 

(31) 

3.1. One dimensional system : Random transverse field Ising chain. 

3.1.1. Gap distribution : finite size analysis. The RTFIC coupled to a ohmic bath was 
treated in detail in Ref. [16]. We just recall here the main results. Since the last spin 
can either be frozen (i.e the last field h is zero) or non-frozen we split Pl^H/Tq) into two 
parts: 

PL{h/To) = ALpLih/To) + (1 - AL)6{h/To) , (32) 

where PLih/Vo) is the restricted distribution of the last fields in the samples that are 
non-frozen and Al is the fraction of these samples. It, or equivalently Pi(log(ro//i)), 
represents the distribution of the smallest excitation energy in the ensemble of non- 
localized spins. At low coupling (small Jq or small a), PLilogiVo/h)) shows indications 
of Griffiths-McCoy singularities characterized by the following scaling behavior for Pl 

PL{\og{To/h)) = V{\og{T,/hU)) , (33) 

where 2; is a dynamical exponent continuously varying with (Jq, etc.). As the 
coupling is increased, z is also increasing and eventually, at some pseudo-critical point, 
Pi:,(log(ro//i)) exhibits a scaling which is characteristic for an IRFP: 

PL{\og{V,/h)) = L-^PiRpp(L-'^ \og{V,/h)) (34) 
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whith ip — 0.32 [16] is a critical exponent characterizing the IRFP. Notice that 
this value of ip is different from V'rtfic = 1/2 computed exactly for the RTFIC 
[3]. The main striking point in the case of ohmic dissipation is that although the 
restricted distribution PL(}og(To/h)) displays Griffith's like behavior like in Eq. (l33ll the 
magnetization becomes finite above a certain length scale L*. This finite magnetization 
is a manifestation of the "frozen" clusters which lead to the concept of rounded quantum 
phase transitions in the presence of dissipation [T4]. Due to these "frozen" clusters, the 
amplitude Al decays exponentially above L*, with Al oc e~^^^* . However, as we pointed 



out in Ref. [16], the interpretation of the finite size analysis (I32|33|34p in presence of 
dissipation has to be done carefully. Indeed in Ref. [E] we suggested that despite 
the presence of these frozen clusters, Griffith's singularities should be observable in the 
susceptibility x{T), above a certain temperature T* oc L*~^ as well as in the specific 
heat Ct,(T). This property can actually be shown (see [Appendix A[ ) on a toy model 
where one considers a RTFIC without dissipation but with a finite fraction p of zero 
transverse fields. 

Here we will use this strong disorder approach to extract thermodynamical 
properties of the full problem described by the Hamiltonian (JS]). 



3.1.2. Susceptibility at finite temperatures. The SDRG successively eliminates degrees 
of freedom with a large excitation energy from the starting Hamiltonian, thereby 
reducing continuously the maximum energy scale of the effective Hamiltonian. If 
continued down to the smallest energy scale the final effective Hamiltonian (consisting 
only of a single but large cluster in an effective transverse field) provides information 
about the ground state of the starting spin chain, the gap, the size, the geometry, etc... of 
the smallest excitation energy. To extract information on thermodynamic properties, at 
low but non- vanishing temperatures, one has to stop the RG procedure at an energy scale 
of the same order of magnitude as the temperature T: clusters (or degrees of freedom) 
that are already frozen at this energy scale will not be active at this temperature and 
behave like classical spins (at this temperature). The thermodynamical properties, 
observables like susceptibility or specific heat, will be determined by the active, i.e. not 
yet frozen clusters. 

It is instructive to have a look at the number and size of frozen clusters as a function 
of the upper cut-off energy, which we identify now with the temperature T. As one can 
see from the left panel of Fig. [2] the number of frozen clusters is zero at high temperatures 
(simply because a < 1) and increases rapidly with decreasing temperature before it 
reaches a maximum and then decreases. The initial increase is due to the formation of 
many small clusters that behave like classical spins at the corresponding temperature 
with moments of the order of 10. The subsequent decrease of the number of clusters 
correlates with an increase in the size of the clusters as can be seen in th right panel 
of Fig. [2] and which is due to the coalescence of small clusters into larger ones at the 
corresponding temperatures. 

With this picture in mind we estimate the zero frequency susceptibility x{T), as 
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Figure 2. Number (left) and size (right) of frozen clusters in the disordered chain 
coupled to an ohmic bath as a function of temperature T for different values of a. 
Shown are data for a single large disorder realization, the size is L — 40000, the 
disorder strength is ho = 1, Jo = 0.025. For this value of Jq the pseudo critical point 
is located at a = 0.2 [16]. 





the sum of two contributions x(^) = Xactive(^) + Xfrozen(^), one arising from the active, 
i.e. non "frozen" spins, Xactivc(^) and one from the "frozen" ones, Xfrozen(^)- In doing 
this, we assume that the interaction between the frozen and the non frozen clusters is 
negligible. Xactive(7') is given by (see also Eq. flA.711): 

Xactive(T) = / -p(e)(l - exp (-/3e)) (35) 
Jo e 

with P = 1/T. To estimate the density of states p(e) using our RG scheme we compute 
the distribution of the amplitudes of fields and bonds which are decimated during the 
renormalization procedure [l]. Having computed p(e), we then perform numerically 
the integration in Eq. (i35l ) to compute XactivciT). In the Griffith's region where the 
restricted distribution scales with L as in Eq. ([33l) . one has p(e) oc e^^+^/^ and thus 
XactivciT) OC T~^^^^^. On the other hand, each (quantum mechanically) frozen spin 
contributes the susceptibility by an amount of 1/T and thus 

/rp\ A/frozen (^) /o^?^ 
Xfrozen(T) = (36) 

where A/frozen (^) dcuotcs the number of frozen spins at temperature T and its finite T 
dependence is computed as explained above. We have computed x(T) using our RG 
scheme for a system of size L = 4096 for different values of a = 0.15,0.17,0.19,0.20 
and 0.22 for Jq = 0.025. In each case, x(T) is averaged over 10"^ different realizations 
of the random couplings and the plots are shown on the left panel of Fig. [3l Let us 
first consider the curves for a < 0.2, where the restricted distribution Pi(log(ro//i)) 
shows a scaling with L as in Eq. ( l33ll [T6| . For low temperatures still above some 
temperature T*, T > T*, one sees in the left panel of Fig. [3] that x{T) is dominated 
by Xactive(T) (X T-i+i/-, for a < 0.2. In this regime of dissipation, one observes 
that the slope of x{T) in a log-log plot depends on a: this is the characteristics of 
Griffith's behavior. At lower temperature T < T*, x{T) is dominated by the 1/T 
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Figure 3. Left: Magnetic susceptibility x(^) for ^ disordered chain coupled to 
an ohmic bath as a function of temperature for different coupling strength a. The 
susceptibility is calculated as x{T) — XactWc{T) + Xfrozen(J") using Eq. (|35l) and Eq. 
([36ll. The size is L = 4096, the disorder strength is ho = 1, Jo = 0.025. For this 
value of Jo the pseudo-critical point is located at a = 0.2 [TQ Right: Crossover 
temperature T* as a function of a extracted from the data on the left panel by the 
condition Xactive(7'*) = Xfrozen(T'*) (see text). 
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behavior oi XfmzeniT) coming from the frozen clusters. Thus Fig. [3]shows that Griffith's 
behavior can indeed be observed above T*. For a > 0.2, the susceptibility behaves like 
~ 1/T in the whole range of temperature. 

In addition, given that we compute separately Xactive(^) in Eq. (l35l) and Xfvozen{T) 
in Eq. ( l36ll our numerical RG procedure allows to estimate T* for which these two 
contributions are equal, XnctiveiT*) = XfmzcniT*). In the right panel of Fig. [3l we show 
a plot of T* as a function of a. One observes in particular that T* shows an inflexion 
point as the pseudo-critical point is crossed such that T* is actually quite small in the 
Griffith's region. 

We now turn to the specific heat Cjj(T) of the spin degrees of freedom. Assuming 
that one can also neglect the interaction between frozen and non frozen clusters one 
immediately obtains that Cy{T) = active (T) given that Cy{rozen{T) = 0. Cv{T) is thus 

C.^'-^ ' ' (37, 



£{T) -£(T = 0) = ^ de Pie) 



exp {-(3e) 



l + exp(-/56) ^^^^ 

where S{T) is the internal energy at temperature T. In the Griffith's region where 
the restricted gap distribution has a finite L scaling as in Eq. (l33ll . one expects 
£{T) — £{T = 0) oc T^+^^^, thus Cy{T) oc T^/^ without any cut-off at some temperature 
T*. In analogy to x{T) we have computed numerically SiT) — S{T = 0) (also averaged 
over 10"^ disordered samples) for different values of a. On Fig. [4] we show a plot of 
T~^{£{T) — £{T = 0)) oc Cy{T) as a function of T for different values of a. One 
observes clearly that the slope decreases as a is increased, i.e. as the critical point is 
reached. We tried to extract an estimate of the dynamical exponent z by fitting the 
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curves in Fig. [4] by T~^{£{T) — £(T = 0)) oc T^/^ at low T as well as by fitting the 
curves in the left panel of Fig. [3] by xiT) oc T^^+^/^ for T > T*. Both estimates for z 
coincide approximately but since the data shown are close to the pseudo-critical point 
(which corresponds here to a = 0.2), it is rather hard to extract properly the dynamical 
exponent given that 1/z becomes quite small, thus one would certainly need smaller 
temperatures to obtain a reliable estimate of z. 

We conclude this paragraph by noting that the data in Fig. [3] and [4] indicate that 
Griffith's behavior of thermodynamical quantities is observable also in the presence of 
dissipation. 
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Figure 4. The scaled energy T~^{£{T) — £{T — 0)), which is proportional to the 
specific heat Cy [T) , for the disordered chain coupled to an ohmic bath as a function of 
temperatur T in a log-log plot. Copmputations are done using Eq. (jM)) . The different 
curves correspond to different values of a. Note that the low temperature behavior is 
characterized by different slopes for different a, corresponding to a varying exponent 
1/z. The system size is L = 4096, the disorder strength /iq — 1, Jq — 0.025. For this 
value of Jo the pseudo critical point is located at a = 0.2 [16] . 



3.2. Disordered ladder. 

Our previous study on Ref. [H] was restricted to the one dimensional case. Here, we 
implement numerically the real space renormalization defined by Eq. (I23tl25p and Eq. 
f l27ll30ll for a disordered ladder coupled to a dissipative bath. When considering a ladder 
(as well as a two-dimensional square lattice) these decimation rules have to be slightly 
modifed to take into account the topology of the system [19]. First, Eq. (l25ll has to 
be modified. In this case, the two spins 2 and 3 are combined to a cluster but when 
we compute the interactions between this cluster and the rest of the chain, one has to 
consider the case in which the two original spins 2 and 3 were actually coupled to the 
same spin i. Although this does not happen in the initial ladder, such a situation may 
occur during later stages of the renormalization. In this case we set the ferromagnetic 
coupling of this spin i with the newly formed cluster to 

<^i,cluster = niax( Jj2, Jiz) ■ (39) 
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The sum of the two bond strengths could also be taken, but does not make a significant 
diff'erence when the probability distribution of the bond strengths is broad. 

The decimation rule on Eq. (I30l ) has also to be modified. This rule says that 
when the spin on site 2 is decimated, effective interactions are generated between the 
neighboring sites of 2. But during renormalization of the ladder there might already be 
bonds Jij present between neighboring sites i and j of site 2. In this case we replace 
Eq. m by 

Jij ~ max(Jij, . (40) 

The topology of the system changes drastically under renormalization. One starts with 
a ladder and the decimations change its structure into a random graph, but this change 
is straightforward to implement numerically. 

In the absence of dissipation a critical point was found for Hq = 1.9, Jq = I [19j. In 
the following we fix Hq = 1.9 and Jq = 0.001 and we vary a. As it was done previously 
for the disordered chain in Ref. \16\ we first focus on the restricted distribution of 
the last fields in the samples that are non frozen (l32l) . For small a, PL^h/To) displays 
Griffith's like behavior as in Eq. ( l33l) . In the left panel of Fig. [5l one plots PLih/Vo) as 
a function oi\og{TQ/hU) with z = 1.7 for different system sizes L = 64,128,256,512 
for a = 0.2. The good data collapse of the curves for different L is in a good agreement 
with Griffith's scaling (i33l ). We observe that the dynamical exponent z increases with 




5 10 15 20 23456789101112 



log(ro/h*L') log(ro/h)/L'^ 

Figure 5. Scaling plot of the restricted distribution pLih/To) of the last field to be 
decimated in the non-frozen samples for the disordered ladder coupled to an ohmic 
bath. Here ho = 1.9, Jq = 0.001. Left: In the Griffiths-region (a = 0.27): pLih/To) 
as a function of \og(To/hL^) for different system sizes L. The best data collapse is 
obtained with the dynamical exponent z = 1.7. Right: At the pseudo-critical point 
{a = 0.27): PLiTo/h) as a function of log (Fo/Zi) /L''' for different system sizes L. 
The best data collapse is obtained with the exponent ijj — 0.27(3). 



increasing a. This is depicted in the left panel of Fig. [6]where one plots again Piih/To) 
as a function of log(ro//ilv^) for different system sizes L but with z = 7.5 for a = 0.262. 
However, despite the fact that the gap distribution displays Griffith's behavior, the 
magnetization is already finite. This can be seen by computing the magnetic moment 
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log(ro/h*L") L 

Figure 6. Left: Scaling plot of the restricted distribution pLih/To) of the last field to 
be decimated in the non-frozen samples for the disordered ladder coupled to an ohmic 
bath. Here ho = 1.9, Jo = 0.001 and a — 0.262, i.e. close to the pseudo-critical point 
at a = 0.27. The best data collapse is obtained by the dynamical exponent z = 7.5. 
Right: Magnetic moment ^(i) as a function of L in the Griffith's region (parameters 
as in the left panel). The linear behavior implies a non- vanishing magnetization rricq 
per spin. 



of the last remaining cluster as a function of the system size L, see the right panel of 
Fig. [H This behavior, which is due to frozen clusters is very similar to the one observed 
for the disordered chain ^TQ\ . 

Finally, one reaches a value of a where z is diverging and one observes a scaling 
characteristic for an infinite randomness fixed point as in Eq. (l34l l. This is shown on 
the right panel of Fig. [5] where we plot L'^Pl (Fo/Zi) as a function of log (To/h)/ L"^ with 
= 0.27(3). 

3.3. Two-dimensional square lattice. 

We have also implemented the decimation rules in two dimensions for a square lattice. 
Here also the topology of the system changes drastically during renormalization. In 
the absence of dissipation a critical point was found for Hq = 5.35, Jo = 1. Here we 
include dissipation, fix Jq = 0.0001 and vary a. At small a one observes Griffith's like 
behavior of the restricted distribution as in Eq. (I33D . In the left panel of Fig. [71 we plot 
PL(\og(TQ/h)) as a function of log (Tq/H) for different system sizes L = 8, 16, 32, 64 and 
a = 0.3. On the right panel, we show that these curves for different fall L on a master 
curve if one plots them as a function of log{TQ/ HL^) with z = 3.1. As we increase the 
value of a, one observes that z is also increasing and eventually we identify a pseudo 
critical point, here for a = 0.37, where the restricted distribution has a scaling form 
characteristic for an IRFP as in Eq. (l34l) with ip = 0.32. This is shown in Fig. [3 
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Figure 7. The restricted distribution PlQi/To) of the last field to be decimated in 
the non-frozen samples for the disordered square lattice coupled to an ohmic bath. 
Parameter values fall into the GrifRths-region: the disorder strength is /lo — 5.35, 
Jo = 0.0001 and the coupling to the bath is a = 0.3. Left: PL{\og{T q / h)) as a 
function of log iTo/h) for different system sizes L. Right: Scaling plot of the data in 
the left panel: PL{\og{T q / h)) as a function oi \ogiT q / hL'^) for different system size L. 
The best data collapse is obtained with z = 3.1. 
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Figure 8. The restricted distribution PL{h/To) of the last field to be decimated in 
the non-frozen samples for the disordered square lattice coupled to an ohmic bath. 
Parameter values correspond to the pseudo-critical point: the disorder strength is 
ho = 5.35, Jo = 0.0001 and the coupling to the bath is a = 0.37. Left: Pl (log(ro//i)) 
as a function of \og{To/h) for different system size L. Right: Scaling plot of the 
data in the left panel: L^PL{log{rQ/h)) as a function oflog(ro/h)/L'^ . The best data 
collapse is obtained with the exponent = 0.32. 



4. Super-ohmic dissipation. 

We have implemented numerically the decimation rules for the super-ohmic bath, which 
corresponds to s > 1. In this case the amplitude A in Eq. (fTS]) and A' in Eq. (l28l) 
which enter the decimation rules are given by 
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A' 



exp 



a2 



1 



1 




(41) 



For s > 1 iterations of the decimation rules (l27ll always converge to a fixed point value 
> given in Eq. ([29l) . Consequently, the spins can not be frozen by the dissipative 
bath. 

We first present results for s = 3, which corresponds to a phonon bath, and one 
fixes the coupling to te bath to a = 0.5 and the strength of the random transverse 
field to Hq = 1.0. All data presented here were obtained by averaging over 10^ different 
realizations of the disordered couplings. 

We first focus on low value of Jq. In Fig. [H one shows a plot of PlITq/H), the 
distribution of the transverse field acting on the last remaining cluster as a function of 
\og{To/hL^) for different system sizes with z = 2.87. The good data collapse of these 
different curves suggests that PL{h/To) exhibits Griffith's behavior: 

Pz.(log(ro//i)) = V{hgiTo/hL^)) . (42) 

Notice that, at variance with the case of ohmic dissipation (l32l l one has here = 1. 
If one increases Jq, z is also increasing and for some critical value of Jq = Jqc here 
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Figure 9. Scaling plot of the probability distribution of the last decimated field for the 
disordered chain with super-ohmic dissipation, here s = 3: PL{To/h) vs. log (Fo/Zil/^) 
(scaling in the Griffiths region) for different system size L. The best data collapse is 
obtained with z = 2.87. The disorder strength is ho = 1.0 and Jq — 0.55 and the 
coupling to the bath is a = 0.5. 



Joc = 0.78 one observes a scaling characteristic for an IRFP 

Pz.(log(ro//i)) = L'^PiKFp{L-'f' log(ro//i)) , (43) 

with ip = 1/2 as in the case without dissipation [3]. This is shown in the left panel of 
Fig. [ini In the absence of dissipation random fields and random bonds play a symmetric 
role in the RTFIC. This is in principle not the case when one includes dissipation in the 
Hamiltonian ([T]). However, this symmetry is restored asymptotically, close to the critical 
point. To show this, we have computed Pl{J /Tq) where J is the last decimated bond. 
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log(ro/h)/ log(ro/J)/L^ 

Figure 10. Left: Scaling plot of the probability distribution of the last decimated field 
for the disordered chain with super-ohmic dissipation, here s = 3: L'^PL{^og{To/h)) 
vs. (log(ro//i))/L'^ (IFRP scaling) for different system size L with ip = 1/2. The 
disorder strength is ho = 1.0 and Jq = 0.55 and the coupling to the bath is a = 0.78. 
Right: The same as in the left panel but for the bond distribution instead of the field 
distribution. Note the similarity of the two distributions. 



In Fig. do] we show a plot of L'^PL(log(ro/ J)) as a function of log(ro//i)) with 
ip = 1/2 for Jo = Jqc- The good data collapse, together with the similarities between 
the plots shown on both panels of Fig. HO] suggest indeed that this symmetry between 
bonds and fields is restored at the critical point. 

To characterize this IRFP, we have also computed the combinations of the products 
of the exponents ipip where ip is another independent exponent associated to this IRFP. 
This can be measured by computing the disorder averaged correlation function C(r) at 
the transition. We compute it by keeping track of the clusters during the decimation 
procedure and compute C(r) = L~^J2iU!i,i+r where Wij = 1 if the sites i and j belong 
to the same cluster, and Wij = otherwise. We have checked that for RTFIC without 
dissipation at the critical point this gives the correct exponent |3j within 5% accuracy. A 
plot of C(r) is shown on Fig. [TT]for different system sizes L = 64, 128, 256, 512 and 1024. 
This plot shows that C(r) oc r~^ with rj = 0.38(1), as in the case without dissipation [3]. 

We have repeated the same procedure for different values of s and found the critical 
value Joc(s). We thus obtain the phase diagram in the plane (Jo,l/s) shown in Fig. 
fT2] where a critical line separates a paramagnetic phase from a ferromagnetic one. 
Along this line, we have found a scaling like in Eq. (I43l) with an exponent ip = 1/2, 
independently of s. One can actually estimate the shape of the critical line in Fig. [T2]by 
assuming that the main effect of dissipation is to reduce the amplitude of the random 
transverse field Hq to given by Eq. (l29ll . If one further assumes that bonds and fields 
play a symmetric role at the critical point (which is fully compatible with our numerical 
results in Fig. [10]), the critical point is then given by the relation Jo = ^o^ the 
RTFIC [3|. This leads to 

Joc = hoexp(—^) , (44) 
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Figure 11. The disorder averaged correlation function C(r) at the critical point for 
the chain coupled to a super-ohmic bath (here s = 3) for different system size L. The 
disorder strength is /iq = 1 and Jq = 0.78 and the coupling to the bath is a = 0.5. The 
error bars are smaller than the size of the symbols. The decay exponent ~ 0.38 gives 
(pij: ~ 0.81 as in the case without dissipation [3J. 
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Figure 12. Phase diagram for a disordered chain coupled to a super-ohmic bath 
(except for s = 1) and ho = l^a = 0.5. The error bars are smaller than the size of the 
symbols. Along the line the critical behavior is governed by an IRFP with -0 = 1/2 
(see Fig. [10] above). The dashed line is the exact expression of the critical coupling 
given in the text l(44|) . 



which is actually in very good agreement with our numerical estimates for the critical 
line in Fig. [121 Using the same arguments, one can also derive analytically the behavior 
of the dynamical exponent when approaching the critical point, this yields 

^ 0^ 777 ' = o log T , (45) 



2(5 ' 2 ^\ Jo 

which we have checked to be in good agreement with numerical results. 

Our results thus suggest that for s > 1, the large scale L > Lg properties of the 
system with super ohmic dissipation behaves at criticality as the dissipationless system. 
On the other hand, one expects that Ls is diverging when s — 1. To estimate its 
behavior close to s = 1 one observes that the typical energy scale at criticality is given 
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by Eg oc /ioexp(— — 1)) ( l44l) . But given that the critical behavior is governed by 
an IRFP, one expects that Lg oc (log-Es)^/'^, with ip = 1/2^ see Fig. [TOl Therefore we 
estimate 

„2 



cx ^ , (46) 



the length above which the system with super ohmic dissipation behaves like the one 
without dissipation. 



5. Conclusion 



In this paper we have developed a real space renormalization, which combines the 
SDRG for strongly disordered quantum magnets with the adiabatic renormalization 
for spin-boson systems, to study disordered, ferromagnetically interacting transverse 
Ising systems coupled to a dissipative bath. In the important case of ohmic dissipation, 
we have first extended our previous study in Ref. |16| to describe thermodynamical 
properties. In particular we have shown that Griffith's Mc-Coy singularities are visible 
in the spin specific heat at all temperatures and in the magnetic susceptibility above 
a (small) temperature T*. For weak dissipation this temperature is extremely small 
and system sizes above which classical behavior in the susceptibility becomes visible are 
extremely large, which represents a major obstacle for numerical studies |21| . 

We have also shown that the disordered ladder as well as the 2d disordered square 
lattice coupled to a ohmic bath displays the same behavior. Using this real space 
renormalization, we also studied the case of super-ohmic dissipation (s > 1). There we 
have found a quantum phase transition described by an IRFP, which is the same as 
the one found without dissipation. Such a scenario is expected to hold also in higher 
dimensions. 

It would be natural to extend this approach to sub-ohmic dissipation (s < 1). 
Unfortunately, it is well known that in that case the adiabatic renormalization fails 
to describe correctly the single spin-boson, which in itself has been the subject of 
recent works [22j. Therefore the problem of an infinite chain (possibly disordered) 
coupled to a sub-ohmic bath remains a challenging problem which surely deserves further 
investigations. 

A final remark concerns the effect of dissipation upon magnetic systems with 
a continuous symmetry instead of the discrete (Ising) case we studied in this work. 
Griffiths-McCoy singularities are much weaker in systems with a continuous symmetry 
[23] . and one would therefore expect that coupling to a dissipative bath would not 
freeze the strongly coupled regions, but enhance their singular behavior. What actually 
happens can elegantly be classified according to whether rare regions including their 
long range interactions in imaginary time due to dissipation are below, at or above their 
upper critical dimension [I5]. A disordered itinerant antiferromagnet, for instance, 
was recently studied with the strong disorder renormalization group and an infinite 
randomness fixed point was found [24], including the accompanying algebraic Griffiths- 



Finite temperature behavior of disordered quantum magnets with dissipation 



21 



McCoy singularities. On the other hand non-intinerant antiferromagnets, involving 
localized magnetic moments, in spatial dimensions larger than 2 like the Heisenberg 
antiferromagnet on the square lattice, will not show pronounced Griffiths-McCoy 
behavior since here the Neel ordered ground state is very robust against disorder [25] 
and no quantum critical point occurs. The effect of dissipation upon strongly disordered 
magnets thus depends crucially on the effect of disorder itself on the system's ground 
state. 
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Appendix A. A toy model for an Ising chain with ohmic dissipation 

To understand qualitatively the full problem described by the Hamiltonian ((31) with 
ohmic dissipation, it is instructive to consider a simpler model where one considers a 
RTFIC without dissipation but with a finite fraction p of zero transverse fields. We 
thus study in detail in this appendix the RTFIC Hamiltonian with k sites having zero 
transverse fields {p = k/L) 

L L 

^ = - E ^^^'^i + E ^i^f and hi, = ... = h,^=Q . (A.l) 

i=l i=l 

First, one immediately sees that the distribution PLih/Vo) shows the same behavior as 
in Eq. (1321 ) with Al ~ e~^l^* and in the small p limit, L* oc p~^ . Besides, the local zero 
frequency susceptibility is 

X.(^ = 0)= f dT{al[T)aim (A.2) 



= f rfrisp{pe«Vfe-^V,n 
JO Z, 

1 ^ — l3Em ^ — f3En 

= ^ E F _ F — 

+ ^ E /?e-^^1(n|ar|m)r, (A.3) 

{n,m {En=Em)} 



where {\n)} is a complete basis of eigenvectors of Ti, (lA.ll) . Their corresponding 
eigenvalues En are such that Eq < Ei < E2 < .... The first term in Eq. ( 1A.3I1 yields at 
zero temperature in the non-degenerate case (all transverse fields positive, finite system 
size L) the known formula 

xrv=o)=2i:Mi|^, ,A.4) 

since then (0|crf |0) = and the last term in Eq. flA.3P vanishes. 
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If one or more transverse fields vanish the Hamiltonian becomes block-diagonal. We 
choose 2;-representation, such that states can be denoted ip = \Si, . . . , Sl), with Si = ±1. 
For convenience we permute the components such that the site ii,...,ik with the k 
vanishing transverse fields /ij^, . . . , hi^, stand to the left: ip = |S'j^, . . . , S'^^,, S'j^, . . . , Sjj^_^,). 
All 2'' blocks are identical up to the diagonal part JiSiSi+i. As a result the two 
blocks belonging to the states with S'j^ = . . . = Si^, (ferromagnetically aligned "frozen" 
spins), have the lowest ground state energy. Obviously 

{^W!Ji') = S^, for p=l,...,k (A.5) 
{i:\at\ij') = for z = l,...,L, (A.6) 
and {Si„...,Si^) ^ {S[^, . . . , S'^J . 
At low temperatures (T — > 0) the main contributions in the sums in (1A.3P comes form 



the terms with En = Eq or Em = Eq, the ground state energy. When k there 
are two ground states lO^), one with S = +1, one with S = —1. For T ^ Z can 
be replaced by 2e-^^o, since Z = 2e-^^o(l + Em e"'^^^™"'^"^). The two ground states 
produce also an extra factor 2 (in addition to the one for the sum over n ^ m, where 
either \n) or |m) can be the ground state): 

4 f,-l3Eo _ p-l3En 9 

= 2E^^ — ^\{nwm\' + Pimm' ■ (a.t) 

-£^« — E- 







The usual argument leading to Xi{^ = 0) cx T^^+^/^ in the Griffiths-McCoy phase of 
the RTFIC with k = involves neglecting the terms n > 1 in the first sum in ( lA.Tll . 
This leads to Xi{^ = 0) ~ (AE)^^, where AE = Ei — Eq is the gap, which follows 
the distribution P{A) ~ A^^+^/^. In the present case this distribution has a cut-off 
at Amin that is exponentially small in L*, the average distance between sites with zero 
transverse fields. Thus one expects the first term of (lA.TI) to display T"^+^/^-behavior 



down to a temperature Tmin = A^J}^ ~ e °^*. The second term is (3 ■ mf, where rrii is 
the local magnetization in (one of) the ground states - and is non-zero due to (1A.5I1 . It 



decays exponentially with the distance x from the nearest frozen site ip-. rrii e "^Z^, thus 
the average over all sites is approximately 



L*/2 Jo 

Thus a T~^-behavior coming from the second term in ( 1A.7I) with amplitude of order 
1/L* competes with a T^^+^/^-behavior with amplitude of order 1 coming from the first 
term in flA.Tp . The latter dominates for temperatures above a cross-over temperature 
Across, which is given by 

T„oss ~ (L*)-' , (A.9) 

which is larger then Tmin (caused by the finite average length L* of the segments) but 
still very small when L* 1 (for instance for L* = 10^ and z = 2 one has Tcross ~ 10^®). 
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